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We report on a preliminary analysis of the diffuse gamma-ray observations of local giant molecular clouds Orion A 
and B with the Large Area Telescope onboard the Fermi Gamma-ray Space Telescope. The gamma-ray emission 
of the clouds is well explained by hadronic and electromagnetic interactions between cosmic rays and nuclei in the 
clouds. Consequently, we obtain the total masses of the Orion A and B clouds to be (80.6 ± 7.5 ± 4.8) X 10 3 Mq 
and (39.5 ± 5.2 ± 2.6) X 10 3 Mq, respectively, for the distance to the clouds of 400 pc and the Galactic CR 
spectrum predicted by GALPROP on the local observations of CRs. The structure of molecular clouds has 
been extensively studied by radio telescopes, especially using the line intensity of CO molecules (Wco) and 
a constant conversion factor from Wco to JV(H2) (= -Xco)- However, this factor is found to be significantly 
different for Orion A and B: 1.76 ± 0.04 ± 0.02 and 1.27 ± 0.06 ± 0.01, respectively. 
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1. INTRODUCTION 

Diffuse emission of > 100 MeV gamma rays in the 
Galaxy is mainly induced by hadronic interactions be- 
tween the Galactic cosmic rays (CRs) and interstellar 
medium (ISM) , via the production of n° particles and 
their subsequent decay into photons. This emission 
can be used to study CRs and the structure of ISM, 
because its emission can be written as the product 
of the CR flux and the density of ISM, and because 
the gamma-ray spectral shape preserves the CR prop- 
erties. Observations of the emission have mainly two 
unique advantages. First, we can study CRs at the lo- 
cation of the gamma-ray emission, while direct mea- 
surements of CRs at the Earth cannot extract their 
positional information which is already lost during 
propagation in the interstellar magnetic fields of the 
Galaxy. Second, ir° emission is not affected by the 
gas condition, such as the temperature of interstellar 
dust, and the CO to H2 ratio which are often used to 
estimate the column density of ISM. 

Since the early stage of gamma-ray astronomy, 
diffuse emission from H I gas and molecular clouds 
(MCs) has been extensively studied to understand 
Galactic CRs and the ISM (e.g. [10]). Among them, 
the Orion A and the Orion B MCs are two of the 
best targets, and have been well studied by pioneer- 
ing gamma-ray telescopes [!, 0, [H, HI ■ This is because 
they are considered to be the archetypes of giant MCs 
located near the Earth, and have been surveyed by 
radio telescopes especially using the line emission of 
CO molecules (e.g. @, HI). From the analysis side, 
their gamma-ray fluxes are strong enough to be distin- 
guished from other diffuse emission components (H I 
gas, inverse Compton, extragalactic diffuse emission), 
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because they located about only ~ 400 pc away from 
the Sun, their total mass is of the order of 10 5 M Q , 
and their Galactic coordinates (£ ~ 210°, b ~ —15°) 
are far from the Galactic plane and the Galactic cen- 
ter. 

In the study of MCs, the biggest difficulty is that 
H2, the main component of MCs, cannot be ob- 
served directly. Hence, CO molecules which are the 
second abundant molecules, have been widely used 
as a tracer of H2 distribution with a conventional 
factor Xco = N(R 2 )/Wco which converts veloc- 
ity integrated CO line intensity, Wco to H 2 col- 
umn density, iV(H2). For example, Xco = (1.8 ± 
0.3) x 10 20 cm _2 (K km s -1 ) -1 was derived by com- 
paring a CO survey with H I and dust observations 
over a Galactic scale 0. This factor has been also 
determined by diffuse gamma-ray observations, e.g . 
(2.3 ± 0.3) xlO 20 cm- 2 (K km s" 1 )" 1 by COS-B\v$, 
and (1.35 ± 0.15)xl0 20 cm" 2 (K km s" 1 )" 1 for the 
Orion region by EGRET 0. However, due to the 
limited angular resolution and photon statistics of the 
preceding telescopes, they could not study Xco for in- 
dividual clouds. Besides, a mystery of "GeV excess" 
which shows a disagreement between the observed dif- 
fuse gamma-ray spectrum and the local CR spectra, 
was reported by the EGRET 0. 

Since the EGRET era, much progress has been 
made in gamma-ray observations and studies on the 
Orion A/B clouds. The most important progress is 
the launch of the Large Area Telescope (LAT) on- 
board the Fermi Gamma-ray Space Telescope (Fermi) 
fllj . Its large effective area (~ 1 m 2 ) and wide en- 
ergy band (20 MeV to > 300 GeV) make it possible 
to study 7T° gamma rays with large photon statistics 
above 1 GeV. In addition, the energy-dependent LAT 
angular resolution is a few times better than that of 
EGRET. As a result, we are able to resolve the struc- 
ture of the Orion clouds on a scale of ~ 1 ° using 
higher-energy photons of better angular resolutions, 
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Figure 1: Gamma-ray (200 MeV to 20 GeV) count map 
of the Orion region binned in 0.25° x 0.25° pixels with 
Hammer-Aitoff projection. The coordinate center and 
the radius of the region are (Z,b) = (211°, -17°) and 20°, 
respectively. In addition to the Galactic plane, three 
molecular clouds, Orion A, Orion B and Monoceros R2 
are visible as bright extended sources. 

and also able to compare the gamma-ray energy spec- 
trum with the predicted one based on the CR spec- 
tra observed at the Earth. ir° gamma-ray emissiv- 
ity of hadronic interactions was modeled in several 
articles using recent results of accelerators and the- 
oretical studies (e.g. [H, [H, [HI), so that the mass 
of the clouds can be calculated backwards from the 
gamma-ray emission and the distance to the clouds 
which was measured accurately from recent MASER 
observations of the Orion nebula in the Orion A 

mm nam. 



2. OBSERVATIONS AND DATA 

The data used in this analysis were obtained in the 
nominal all-sky survey mode of the LAT between 4 
August 2008 and 15 August 2009. Among all the LAT 
events, we selected ones classified as P6_V3-DIFFUSE 
class . The reconstructed energy range and zenith 
angle were limited to 200 MeV - 20 GeV, and < 105°, 
respectively. Gamma rays in a circular region of ra- 
dius 20° centered at (£,&) = (211°, -17°) were then 
selected for later analyses. 

All events in the selected data are binned in 160 x 
160 equal area pixels with the Hammer-Aitoff projec- 
tion as shown in Fig. [T] The Orion A and B are easily 
recognized by their count excess from the surrounding 
diffuse emission and point sources. 

The reconstructed energies of the events are binned 
in a logarithmic series of 16 between 200 MeV and 
20 GeV. The LAT exposure of each energy bin was 
calculated from the LAT pointing history and instru- 



mental response function (IRF) using ScienceTools 
v9rl5p4. 



3. ANALYSIS 

3.1. EXTRACTION OF GAMMA-RAY 
EMISSION FROM ORION A AND B 

To study the gamma-ray emission associated with 
the Orion MCs (molecular gas), other emission com- 
ponents must be subtracted from Fig. Q] The emis- 
sion in the region consists of several components. The 
dominant ones are diffuse 7r° gamma-ray and elec- 
tron bremsstrahlung emission induced by interactions 
between CRs and HI/H2 gas. In addition, inverse 
Compton (IC) scattering of CR electrons off interstel- 
lar radiation fields exits in the region. Extragalac- 
tic diffuse emission, and residual instrumental back- 
ground are classified as diffuse emission. We refer to 
the latter two components as "isotropic component" . 
The LAT detection of some point sources in the region 
has been reported. However, there is no intense point 
source overlapping the Orion MCs [l9| . 

In order to extract the Orion A and B, we mod- 
eled other emission components using a computer pro- 
gram called GALPROP which calculates CR flux and 
diffuse gamma-ray emission in the Galaxy [2(| l2lj . 
It uses the data of H I and CO surveys to calcu- 
late the target mass of CR interaction |{J [2a| 1 ■ Its 
calculation is known to be consistent with the LAT 
gamma-ray observations in a galactic scale except for 
the Galactic plane, when reasonable input parame- 
ters based on many observations and experiments are 
given [23|, [24|, HE|. Each set of input parameters is 
referred to as a GALDEF file. In this analysis, we 
use GALDEF 54_77Xvarh7S in which CR proton and 
electron fluxes are scaled by factors of 1.15 and 1.75, 
respectively, to reproduce the LAT observations. 

Before we subtracted the H I contribution in the re- 
gion, its normalization was estimated in nearby two 
10° radius regions near the Orion MCs where no 
known large MCs exists. The regions are centered 
at (£, b) = (230°, -16°) and (210°, -32°), and the ob- 
tained fluxes are 1.15 and 0.98 times those predicted 
by GALPROP, respectively. Thus we multiply 1.07 
to the GALPROP prediction for gamma-ray intensity 
from H I gas, and add 8% to the systematic uncer- 
tainty of H I subtraction process. The scaled H I emis- 
sion was convolved with the LAT IRF, and subtracted 
from the pixclizcd data in each energy bin. In addi- 
tion to this, IC emission was also subtracted. Since 



1 We assumed a constant spin temperature of 125 K for the 
H I map. 
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(a) Extracted gamma-ray intensity (b) Modeled W^co - based intensity 

Figure 2: (a) Same as Fig. [T] but shown after subtracting the other diffuse components. The unit is integrated photon 
intensity between 200 MeV and 20 GeV. Solid lines show the definitions of the regions of Orion A and B. Rectangles 
with dashed lines show the regions which are used to estimate the isotropic component. A straight dot-dashed line at 
I — 212° is the boundary of two separated regions in Orion A. (b) A modeled gamma-ray intensity map simulated from 
a constant Xco = 1.5xl0 20 cm _2 (K km s -1 ) -1 , Wco map 0, the LAT response, the emissivity of 7r° gamma and 
electron bremsstrahlung, and the CR flux predicted by GALPROP. (These figures are still PRELIMINARY) 



the IC contribution is smaller than other emission, its 
uncertainty does not affect our results. 

Finally, the isotropic component was estimated 
from surrounding regions (background regions) 
around the Orion MCs, then subtracted from the 
region-of-interest. Fig. [2a| shows the gamma-ray in- 
tensity map after the subtraction method described 
above . It was divided by the LAT exposure, and re- 
binncd in 1° x 1° pixels. We defined the boundary of 
Orion A and B with solid lines in Fig. [5aJ 



3.2. ENERGY SPECTRA AND TOTAL 
MASSES OF THE CLOUDS 

The energy spectra associated with the Orion A 
and B clouds are shown in Fig. [3] All photons in 
the boundary regions in Fig. [5a] are collected for each 
cloud. They arc fitted with 7r° and bremsstrahlung 
components with two free normalization parameters. 
Here, we calculated 7r° gamma-ray emissivity using a 
recent parameterized model 12j. Input CR spectra 
at the Orion region 2 predicted by GALPROP are as- 



sumed, which are ~ 8% smaller than the observed CR 
flux at the Earth [2(| H3, [H|. Gamma-ray inclusive 
cross sections for proton-He, alpha-H, alpha-He, and 
heavier CR metals are scaled using the method de- 
scribed in [l4|, HH- A factor 1.02 as a contribution 
from ISM metal was finally multiplied to the calcu- 
lated 7T° emissivity [l3l ]. The bremsstrahlung spec- 
trum was calculated using GALPROP. 

The obtained fit results show good agreements with 
the 7T° spectrum. This means that gamma-ray emis- 
sion from MCs is dominated by known physical pro- 
cesses, i.e. interactions between CRs and the gas. 
Therefore, we are able to calculate the total mass of 
the clouds. Assumin g th e distance to the Orion clouds 
to be 400 pc [lE EE ll3, E3I , the mass of the Orion A 
and B are estimated to be (80.6±7.5( S ^ a ^-,±4.8^jj jj) x 

1O 3 M and (39.5 ± 5.2 (stat) ± 2.6 (H T) ) x 10 3 M Q , re- 
spectivcly. The second terms indicate statistical er- 
rors of the fit, and the third ones are systematic er- 
rors introduced by the uncertainty of H I subtraction 
explained in subsection 13. II 



3.3. CORRELATION BETWEEN 
GAMMA-RAY AND CO INTENSITY 



2 Corresponding to a cylindrical Galactic location (R 
8.5 kpc,Z = -0.14 kpc). 



Assumption of a constant Xqq in small scales (~ 
10 pc to ~ 100 pc) or Galactic scale has been widely 
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Figure 3: Energy spectra of the Orion A (black) and B (red) fitted with modeled spectra of it gamma (dashed) and 
electron bremsstrahlung (dotted). Statistical errors are shown with bars, and the systematic errors of the LAT response 
are shown with shaded area. Polygons of solid lines show the systematic uncertainties of H I subtraction process. (This 
figure is still PRELIMINARY) 



used in studies of MCs (e.g. @, 0), while gradient of 
Xco is also discussed [30,|31(. However, there is room 
to reconsider this simple assumption that Woo can be 
used to trace the structure of MCs [H, [33| . 

Utilizing the good angular resolution and large pho- 
ton statistics of the LAT, the correlation between 
gamma-ray intensity and a Wco map can be stud- 
ied in a scale of ~ 1°. Fig. [2b] shows a mod- 
eled gamma-ray intensity map based on a CO survey 
and GALPROP calculation with a constant Xco of 
1.5 xlO 20 cm" 2 (K km s -1 )" 1 , where the LAT IRF is 
convolved. 

Fig. 2] shows the pixel- by-pixel correlation between 
the observed gamma-ray intensity (Fig. l2al x) and the 
Wco model (Fig. I2b[ y). The best fit results by a lin- 
ear function (y ~ Po + Pix) are also shown. If the 
gamma-ray intensity can be presented by the product 
of a constant CR flux and a constant A"cOj the slopes 
of the best-fit functions become constant. However, 
we found that the obtained slopes are significantly dif- 
ferent for the Orion A and B clouds, while these clouds 
are thought to be in a common environment since their 
birth. In Fig. [4^,, two additional linear functions are 
shown: one is fitted with the larger longitude region 
in Orion A {£ > 212°), and the other is for the smaller 
longitude region (212° > t) (see Fig. 121)). 

Since we modeled the vertical values using a con- 
stant X co of 1.5 xlO 20 cm~ 2 (K km s -1 ) -1 , the cor- 
responding Xqo °f Orion A and B are (1.76 ± 



0.04 (stat) ±0.02 (H T) ) x 10 20 and (1.27 ± 0.06 (stat) ± 
0.01,jjj0 x 10 20 , respectively. 

4. DISCUSSION 

We obtained the energy spectra of the Orion A and 
B, and showed that they can be explained by the CR 
interactions with the nuclei in the gas. Thus, a lin- 
ear correlation between gamma-ray intensity and the 
column density of the clouds are expected, because 
gamma-ray emission is not affected by environmen- 
tal condition of the gas. In fact, the correlation be- 
tween gamma-ray intensity and the VFccrbased mod- 
eled map holds a linearity for roughly one decade with 
only a few 10% deviations. This implies that cosmic 
rays of energy above ~ 1 GeV can penetrate dense 
cores of molecular clouds. However, the correlation 
slopes in the Orion region was found to be not con- 
stant as shown in Fig. @] 

There are some possible interpretations of the dif- 
ferent correlation slopes. In the analysis, we assumed 
a constant CR flux in the Orion region. However, if it 
is significantly different in the three separated regions 
in Fig. I2al the gamma-ray intensity also varies accord- 
ing to the CR flux variation which might be caused 
by the strong magnetic field in molecular clouds. 

On the contrary, if the CR flux is almost constant 
in the region, we need to consider nonuniformity of 
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Figure 4: Correlations between the measured gamma-ray intensity (x) and a model prediction by Wco {y) assuming 
Xco = 1.5xl0 20 cm _2 (K km s -1 ) -1 . Data points are fitted with a linear function y = po +pix. (a) Orion A. Circles 
correspond to the data in pixels of I > 212°. (b) Orion B. (This figures are still PRELIMINARY) 



X C o in the region. While the CO line ( J = 1 — 0) 
is the de fact standard of mass tracers of molecular 
clouds, the CO to H2 ratio can be varied by the con- 
dition of each cloud. By comparing Wco maps with 
dust observations or gamma-ray observations [321 l33j , 
it is known that there exists gas which is not traced by 
CO or H I observations. The nonuniformity of Xqo 
in the Orion region can also be explained by this. In 
fact, IR emission from interstellar dust and visual ex- 
tinction by dust are stronger in the left half of Orion 
A than that expected from CO observations [HI, Hf|. 
Therefore, if there exist H 2 molecules not fully traced 
by Wco, but traced by gamma-ray and dust, the Xco 
variation shown in Fig. [4] can be explained. In addi- 
tion to H 2 , it is possible that a part of H I gas was not 
subtracted adequately, because its spin temperature 
is relatively low in cold molecular clouds compared 
to surrounding diffuse H I region. While we assumed 
a constant spin temperature of 125 K in the region, 
optically thick H I gas whose temperature is low may 
not have been corrected properly and contribute to 
the gamma-ray emission, especially in the left part of 
Orion A. 

While the gamma-ray emission from H I and H2 
cannot be distinguished from each other, our observa- 
tions revealed non-linear relation between gamma-ray 
and Wco maps, which will enable us to understand 
the nature of molecular clouds in detail, in addition to 
the Galactic CRs. More photon statistics in future is 
expected to show the smaller structure of the clouds. 
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